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Abstract: We explore the small-temperature regime in the deconfined phase of massive 
fundamental matter at finite baryon number density coupled to the 3+1 dimensional A/" = 4 
SYM theory. In this setting, we can demonstrate a new type of non-trivial temperature- 
independent scaling solutions for the probe brane embeddings. 

Focusing mostly on matter supported in 2+1 dimensions, the thermodynamics indicate that 
there is a quantum liquid with interesting density-dependent low-temperature physics. We 
also comment about 3+1 and 1+1 dimensional systems, where we further find for example 
a new thermodynamic instability. 
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1. Introduction 

In recent years, the ADS/CFT correspondence [1, 2, 3] has become a powerful tool to 
study various properties of the strong coupling limit of conformal field theories, with ap- 
plications to QCD and also more recently to some aspects of condensed matter physics. 
In principle, one has to differentiate between top-down setups that are constructed within 
string theory and imply consistency - and bottom-up setups, where the gravitational duals 
are constructed from a phenomenological point of view. In this paper, we use the former 
approach as we would like to explore what happens to a particular consistent theory. 

In condensed matter applications, there has been particular interest in 2+1 dimensional 
systems that can be typically constructed using M2 branes [4], or as a defect in a 3+1 
dimensional background using D3-D5 [5, 6] and D3-D7 [5, 6, 7] intersections - and also 
as bottom-up setups in various contexts such as superconductivity. As our world is 3+1 
dimensional, the defect setup may be more realistic, even though there are some problems 
with the consistency of the D3-D7 systems [6, 5]. Recently, there has also arisen significant 
interest in fermi- liquid-like aspects at low temperature compared to the density, for example 
in a 3+1 dimensional D3-D7 setup [8] or a 1+1 dimensional D3-D3 configuration [9]. This 
low-temperature limit will also be in the focus of this paper. 

A common way how to introduce fundamental matter in ADS/CFT is in probe brane 
configurations, where one considers an AdS blackhole background. For example one con- 
siders the well-known AdS^ x (black hole) solution (above the deconfinement phase 
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transition) from the decoupling limit of a stack of A'^c ^ 1 D3 branes, that is dual to a 
(thermal) SU{Nc) A/" = 4 supersymmetric Yang- Mills theory [2]. Then one inserts A'^^ <C Nc 
intersecting Dp branes, giving Nj families of (charged) fields in the fundamental represen- 
tation of the SU{Nc), living along the directions of the intersection [10, 11] - such as the 
above-mentioned D3-Dp intersections. 

In these setups, the mass Mg of the fundamental matter is usually given by a scalar 
describing the size of the compact sphere of the probe brane geometry inside the [10]. 
Typically large masses or small temperatures correspond to a vanishing size of the compact 
sphere away from asymptotic infinity, i.e. to a "narrow funnel" that depends for example 
on the baryon number density p. However it was noted in [12] that in d + 1 dimensional 
probe brane setups there is a coincidence in the scalings at finite p and M^: The relevant 
parameters are dimensionless with appropriate powers of the temperature and at small 
temperatures the dimensionless mass scales as ^ ~ a i^^^^^, where the proportionality 
constant a depends on the "size of the funnel" . Hence, the exponents of T cancel and the 
size of the compact sphere remains finite and independent of temperature. This can be 
visualized in the usual picture by the extra tension or "stiffness" of the probe branes from 
the increasing as T decreases, which will be illustrated in figure 1. 

This scaling coincidence suggests that there may be some interesting low-temperature 
physics in this limit, and hence we will explore the details of the embeddings and their 
thermodynamics and related physical properties. We will mostly consider the 2+1 dimen- 
sional defect setup, but eventually we will also briefiy look at the results in 3+1 and 1+1 
dimensions, and point out some important differences. 

This paper is organized as follows: First, we review the string theory setup and the 
necessary ADS/CFT dictionary of the D3-D5 defect in section 2. Then, we will demonstrate 
the scaling solution for the probe-brane embedding that is central to this paper in section 
3 where we also comment on the relevant thermodynamic variables to consider in the 
massive case. In section 4, we first obtain the leading terms of the free energy, then 
discuss the thermodynamics/statistical mechanics in 2+1 dimensions and finally study the 
subleading terms that give rise to the heat capacity. Throughout this section, we will 
see how the system interpolates between a mass-dominated limit and a density-dominated 
limit. Finally, we will briefly discuss the generalized 3+1 and 1+1 dimensional systems in 
section 5; and then we will discuss the results and conclude in section 6. 

2. Setup 

We start with the supergravity background of a planar black hole in AdS^, 

ds' = ^ (-(1 - r'o/r')dt^ + d^) + ^ (^^^ + r^d^l^ , = ~ (2.1) 

This corresponds to the decoupling limit of Nc black D3-branes dual to A/" = 4 SU{Nc) 
super- Yang-Mills theory at flnite temperature T, living along the flat directions of the AdS, 
[3] . The temperature T is given by the Hawking temperature T = and the Yang-Mills 
coupling by = iirgg. Since the curvature L is given in terms of the string coupling Qs 
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and string length Is as = At: g^Nc Ig, the 't Hooft couphng A = Qym^c can be written as 
\ = j^. Hence the "supergravity hmit" L ^ Ig in which the type IIB supergravity action 
and the solution (2.1) are valid corresponds to strong coupling A ^ 1. 

In practice, however, we will use coordinates that are made dimensionless with factors 
of ^ = denoted by (•), such as i := ^ = t (vrT), and the inverse dimensionless radius 
n := ^, giving us 

ds^ = [-{I - u')dP + dil + + dnl^ . (2.2) 

Because in this setup all the fields transform in the adjoint representation of the 
SU{Nc) its use in QCD or condensed matter physics is very much limited as one would 
like to consider also matter that is charged under this symmetry, i.e. that transforms in 
the fundamental representation. To introduce the fundamental matter one then creates 
an intersection of "probe" Dp branes with the D3 branes, such that at the string theory 
side there are fields at the massless level of field theory at the intersection. From the point 
of view of the probe branes, they correspond to endpoints of D3-Dp strings, and in the 
gravity side they correspond to fundamental fields in the (defect) field theory. 

Here, we use the well- known D3-D5 defect setup (see e.g. [13, 14, 15]): 

12 3 

t X y z 

background : D3 x x x x 
probe : D5 x x x 

The dual field theory is now the SYM gauge theory coupled to Nf fundamental hyper- 
multiplets, which are confined to a (2+l)-dimensional defect. This construction is still 
supersymmetric, but the supersymmetry has been reduced from A/" = 4 to A/" = 2 by the 
introduction of the defect. In the limit Nf ^ Nc, the D5-branes may be treated as probes 
in the supergravity background, i.e. we may ignore their gravitational back-reaction. 

Here we want to turn on only the overall U{1) factor of the world-volume gauge field, 
so the the probe branes are governed by the DBI action 

S = -nNf [ ^-det{P[G] + 2TTpgF) , (2.4) 

evaluated in the D3 background (2.1). We assume also translational invariance along the 
flat directions and rotational invariance on the sphere. Hence, the pullback in (2.4) gives 
us one scalar field corresponding to the position in the z direction, which was extensively 
studied in [6, 5], and another scalar which describes the size of the compact sphere and 
corresponds to turning on the mass of the fundamental matter, studied in [12, 5] and more 
extensively in the similar D3-D7 system in [16, 17, 18, 19, 20]. Parametrizing the as 
dQl = dO"^ + sin^ 6 d^ll + cos^ 9 dQ^ and putting the brane on the first S'^ of the S^, the 
induced metric on the probe branes is given by 

ds^ = 5 (-(1 - u^)dP + dxl + + du^ + u\l- nuf)dni) , 

(2.5) 



5 6 7 8 9 



X X 



(2.3) 
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where we defined the scalar as ^{u) = sin 6{u). 

We also choose to turn on the world- volume U{1) gauge field which is dual to the U{1) 
current operator of the U{Nf) that gives rise to the "flavor symmetry". In particular, the 
flux of the electric field 

F = duAt{u)duAdt , (2.6) 
corresponds in the field theory side to the baryon number density (see e.g. [21]) 

P = (J') = -TT^TJ- = ^TrNf^^ lira duAt{u) . (2.7) 
In this parametrization, the DBI action becomes 
S = ATTl'^nNf j daV-detP[G]4Vl " (27r/2a„^t('u))2p[G]**P[G]«" (1 - ^2) (2.8) 

where the integral is taken over the AdS^ part of the world-volume. Using this action, it 
is straightforward to obtain the solution for the gauge field 

= VAT P^^-^^(-) + u^WWW 
^ ' Vl - ^{uY^p^u^ + (1 - ^(n)2)2 

where p := j^jj^- The equation of motion for ^(n) becomes 

2(1 - ^2)3 + ^2(1 _ ^4')(^2y4 ^2)2)^/2 



4(1-^-2) /(l_^2)(l_^2 + (^^2_ y6)^/2) (l + ^2^4 + ^2(^2_ 2)) 



W2 Y (l_^2)(l_^2 + (^2_^6)^/2)^ ' ^^-^i 

and on the horizon u = 1, the equations reduce to 

lim du"^ = 7^:3;r^ , where *o = lim I- . (2.10) 



u-,1- " 2p2 + (l_^2)2 



The asymptotic solution at u — )• is 

^ ~ mu + cu^ , (2-11) 

where m and c are dimensionless free parameters that are determined by the value on the 
horizon and (2.10) that acts as a second boundary condition on the horizon. Following 
arguments of the T-dual (3 -|- 1) dimensional D3-D7 setup [16, 17, 18, 19, 20], the quark 
mass Mq and dual condensate C are given by 



Mg = ^ = ^Tm and C = T^NfM^c . (2.12) 



One can understand this identification of the mass from the separation between the D3 
and D5 branes in flat space, such that Mq is the mass of a stretched D3-D5 string and the 
condensate is just the thermodynamic dual of the mass. 
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"b. 5. 10. 15. 20. 25. (1/u) cos(e) 

Figure 1: The full "physical embedding" of the probe brane, represented by the "radial" direction 
of the bulk spacetime ^ cos 9 and the "separation" ^ sin 6 for fixed p = 1 (at fixed temperature 
and varying mass m, black) and fixed p = ^ = 1 (at fixed mass and varying temperature T = i, 
gray). The corresponding curves of both cases should meet at asymptotic infinity, ^ cos 9 — >■ oo 



At vanishing density, there is a critical temperature-mass ratio below which the probe 
branes do not extend down to the horizon [22, f6]. At finite densities, and unless one turns 
on the scalar in the z direction considered in [5, 12], the brane always extends down to the 
horizon even though a phase transition may still be observed at small densities [12]. In the 
limit considered in this paper, however, this phase transition is of no concern. 

In general, equation (2.9) has no analytical solution and (2.10) implies that one has 
to start integrating the equation from the horizon to obtain the mass and the related con- 
densate at the asymptotic boundary, rather than setting either of them first. We show the 
full numerical solutions for ^(u) in figure 1, where we illustrate the "physical" embedding 
cos 6 , ^ sin 9^ of the probe brane for both fixed p/T"^ and fixed p/M^ to support the 
physical interpretation outlined already in the introduction. 

In the following section, however, we will demonstrate a particular case in which the 
equation of motion for ^ simplifies dramatically. 



3. Taking the limit 

In this section, we obtain the scaling solution in a "blow-up region" near asymptotic infinity 
— 7- at small temperatures and finite density and mass, i.e. p^l, m^l. 

Before we do so, we need to first consider an approximate solution at finite values of u 
in order to match the asymptotic solution to the black hole horizon. Hence, let us consider 
a linear expansion around an approximately finite size funnel, ^{u) = ^'o + V'(^)) which 
reduces (2.9) to 

2*o(l - *o) = ((1 - u'')dui^iu)) + 0{ij^) . (3.1) 

Enforcing the boundary condition (2.10) and ip = at the horizon gives the solution 

-^'o(l-^'^)/ 4 4 il-u 4, (1 -Fn2)(l + n)2\ . 
^(n) = (l - - 2u^ tan- ^^^uHn ^ j + o(,-) , 

(3.2) 
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which is negative over u g]0, 1[ and becomes finite at small u ~ 0{l/^/~p) as it should. 

Now we wish to find some overlap of this solution with a solution in the blow-up region 
at small u that we can match around u 1. To do so, we define the "blow-up" 

coordinate ^ := u^/p, which maps the horizon to large ^ = ^/p ~ T^^ S> 1, and the scale 
over which most of the change in ^ takes place to ^ g]0, 0(1)]. Then, we consider all of 
^ G M^" and expand for large p. The leading term (up to order /0~^) is written in the most 
compact form as 

= 8j (\/l - "PKjV?" + (1 - (^)) P-3) 

where the erratic solution ^' = m ^ has to be ignored. It is straightforward to show that the 
expansion for large ^ is indeed consistent with = ^'o — ~^^~^Jo^^^ i-^- with the small-u 
limit of (3.2). In practice for the numerics, we have to choose as a coordinate and 
implement this matching as a boundary condition as — )■ 0. 

Because we want to consider fixed mass, but vary the temperature and density, the 
suitable physical quantities to consider at finite mass are 

f = - = , (3.4) 

m 2Mq ' ^ ^ 

P = = ^T^T^ and (3.5) 
- ^ A = ^ C a) 

where the numerical coefficients are chosen to give us parameters that are straightforwardly 
related to the embedding. Using (2.7), (2.11) and the definition of ^, they are now related 
to the scalar by 

1 



p = (S^^lg-^o)"' , c = p-dl^\^^o and T = y/p/p . (3.7) 



2 

It turns out that there seems to be an exact relation 5|^'|g_s.o = ~2^0) such that c = — p^o- 
For small ^'o, we can expand the equations of motion to linear order, 

= 2^- 2^9^^ + e^(l + ^^)9|^' . (3.8) 

This has an analytic solution in terms of the elliptic integral of the first kind, J-'{(j)\k) = 
fn I '^^ ^ and the complete elliptic integral }C{k) = Ti'^\k) 

yl— fc^sin ((9 

m = q>QS^(e'^'^F{ishi\i~^{e'^'^i)i) + /C(l/V2)) , (3.9) 



which is itself not particularly interesting because it just interpolates between the linear 
behavior at small ^ and the constant ^'o at large ^; and since it is the linear expansion in 
it tells us that at large densities the relevant scales in ^ do not depend on the density. 
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Figure 2: Left: "foCp) with the approximations at small and large p. Right: for various 

values of p. 



The profile can be seen m the curves with large p on the right in fig. 2. The important 
point is that because VA; : F{x\k) ~ 2: + 0{x^), we can straightforwardly read off the first 
and second derivatives at ^ — t- to find that 

^'o ^ ^ and c ~ -^'oP ~ — ^ r- ■ (3.10) 

In the regime of ^ — 1 ^ 1, however, one cannot find such a simple expression. By 
studying the equations of motion, one can only argue that the transition from the large-.^ 
solution to the asymptotic solution occurs at ^ ~ 0(1 — ^'o)*-^^'*\ such that 1 — \I' oc p^. 
The corresponding value of the condensate is c = —p + 0{p^) for small values of p. 

The full numerical dependence ^q{p) and these limiting cases are shown in figure 2, 
where we also show the corresponding embeddings ^'p(0- We further illustrate^ the "size" 
of the sphere, cosO = ^/Y^^W^ in figure 3. 

4. Thermodynamics 

In order to study the thermodynamic quantities, we follow the standard procedure of 
obtaining the free energy at fixed quark mass as a function of density and temperature, F = 
F^Mqip^T), and then computing its derivatives. The extrinsic thermodynamic quantities of 
the defect should be considered as a contribution to the overall 3+1 dimensional system. In 
[12], it was argued, however, that they can be discussed independently, since the asymptotic 
properties of the bulk do not depend on the details of the defect. As usual, we consider all 
extrinsic quantities densitized and will comment on obtaining the pressure separately. 

4.1 Free Energy 

In [12], we computed the free energy density for this system via the standard procedure 
from the Euclidean action, Jg, [16] using 

F = Tie , ie = Ce + hdy + M , (4-1) 



^The color scheme is an artist's impression to illustrate the similarity with the flowers of Angel's Trumpets 
(Brugmansia). 
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Figure 3: Illustration of the size of the temperature- independent size cos e = ^/l - *2 of the 
compact S"^ (only shown here) in the asymptotic blow-up region as a function of the inverse 
radial coordinate ^, for — 0.95. The color scheme is an artist's impression to illustrate the 
similarity with the flowers of Angel's Trumpets {Brugmansia) . 



where the boundary terms 

hdy. = -IV^+I'^^VI (4-2) 

are dictated to us by consistency [23]. The Legendre transformation needs to be done in 
order to obtain p as a thermodynamic variable, rather than the chemical potential (see 
e.g. [19, 12] for a more detailed discussion). At large temperatures, i.e. m,p <^ 1, this 
evaluates to 

" " ^ -It^ , (4.3) 



with the corresponding entropy S = T'^ and heat capacity cy = . Note that the unusual 
factor A arises because of the definition of T. 

To compute the value of the free energy in our case, we can again use ^ = U\fp and 
expand for large p to obtain 




3^2 



6^^ 



bdy. 



This has to be evaluated over ^ g]0, y/p], and we will try to split it into an integral over 
^ G]0, oo[ and ^ G]\/p,oo[, for the latter of which we can consider a simple expansion of the 
integrand. At large ^, the integral becomes just f d^{^~'^ + 0{£,~^)). To see that this is 
consistent at ,^ ~ 0{^/~p) with the full expression for the integral found in [5], we substitute 
the approximate the solution (3.2) into the full expression and expand for large p - giving 
us / du{^ + 0{p~^)). In contrast to the procedure in [12] and the "main" part of the 
integral ^ €]0,oo[, we did not add the boundary term to this part of the integral. This is 
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because it vanishes in the fuh solution on the horizon and in our approximate solution as 
— )• oo, but not sufficiently fast near ^ ~ ^/p. 

Putting both parts of the integral together, we find the result 

F = p''H{-p) -f-p + 0{f^/VP) , (4.5) 

where the term Tp comes from the integral over ^ &]y/p, oo[ and the integral 



Hp) = lim ^ + (1 - *2)2^i + i!^ - 2 + _ 2^^^/ 



(4.6) 



is a function of p only. An important observation is that the subleading term is of higher 
order in temperature than the behavior 4.6 in the high-temperature regime. 

At ^' = 0, i.e. at p oo, the integral simplifies to tig ^'^^"^ = |/C(l/\/2). We 
can also expand the integral to order ^'^ using the approximate solution for 'f, equation 
(3.9), but the resulting integral can only be evaluated numerically and yields 

m = \k{1/V2) + ^ + Oip-') ■ (4.7) 
3 p 

In the opposite limit at p <C 1, we approximate the solution by two integrals, below 
and above the "singular" point around ^/p, which is the "kink" in the profiles shown in 
figure 2. The integral from to ^/p with * ~ S,/\fp gives a contribution +0(1) and 
the integral from 1/y/p to oo gives ^^3/2 ~^ ^ ^'-^ leading order in 1 — ^q. Hence the leading 
behavior at small p is <I> ~ ^ or F ~ p{l — T). The next order in ^{p) is numerically 
approximately 1.95/3^/^, i.e. 1.95p^ in F. 

We show the full solution for $(/o) together with the approximations in figure 4. 

Reinstating the dimensionful coefficients, the free energy density becomes 




F = ^NrP J^^(p) - T 



N,rhoM,-^N,rhoT : ^ ^ 

(4.8) 
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4.2 Thermodynamic quantities 



The entropy and internal energy densities are given by S 
which evaluates in our case to 

S = ^Ncp + Oity^) and 



III andU = F + TS, 



U 




NrVho Mr, 



Correspondingly the specific heat capacity cy = ^ I . is 



dU \ 
dT\p 



+ 0{f^/^) . 



We can also obtain the chemical potential p 



d_F 
dp 



and the "inverse density" of states e ^ 



^NcT 



(2M, 



d p2 



^NcT 



(4.9) 



(4.10) 



(4.11) 



(4.12) 




+ 0.73 



NfM^ 



« m 
■ ^ » (^)^ 



(4.13) 



where one has to keep in mind that P = ^ and the subleading term in these expressions 
is of order T^. 

The density scaling of the chemical potential in the large-/? limit is in principle char- 
acteristic for a 2 dimensional Fermi liquid. Interpreting this chemical potential however as 
a Fermi-sea -type chemical potential may be misleading, as the results of [5] suggest that 
the finite density gives rise to some finite length scale and to a quasiparticle spectrum. 
Furthermore, for a Fermi sea one would expect a heat capacity proportional to the tem- 
perature. Hence, it should be rather interpreted as a ground state energy that depends on 
the density. 

Somewhat similarly, the entropy and the —^/XNcT/2 term in the chemical potential 
suggest a ground state degeneracy proportional to the total baryon number, or \/A/2 times 
the quark number. It may be an interesting exercise to write down a distribution and 
density of states that reproduces this. 

The interpretation of the system as being in a ground state is also consistent with the 
large-mass limit. There, the energy of the system is approximately the baryon number 
times the baryon mass scale NcMg, or precisely the quark number times the quark mass. 
The cross-over between the high and low density limits occurrs as the induced length scale 



IN 



j- is of the order of the quark mass. 
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To obtain the pressure, we have to keep in mind that F is densitized. Using dyG{p)\N = 
'rdpG{p) for some function G that only depends on /), we can obtain P = —F + pdpF: 




+ 0.59^ 

i/C(^)iV.p3/2, 



• Nt 



■ N, 



> 



(4.14) 

From the pressure, we can then obtain the isothermal and adiabatic bulk moduli Kt and 
Ks- Since the pressure is an intrinsic quantity, we can compute the volume derivative as 



Kt 



-V 



dP 

dv 



T,N 



pd-p {-F + pd-pF) = p^d^pF: 



Kt 




+ 1.76 



lr(J-\N £.3/2 [A 



(4.15) 



and because the full, extrinsic, entropy is just given by the baryon number, we have Kt = 
Ks- Overall, these results suggest that at small densities compared to the mass, we are 
dealing with some type of pressureless gas of quarks, where the pressure scales proportional 
to and p^, compared to pT for an ideal classical gas. As the density increases above 
the mass scale, the zero-temperature pressure and bulk modulus become large. This is 
consistent with the density dependence of the ground state energy and with a length scale 
induced by the finite density. 

We can also obtain the speed of sound. For the small-density limit this is just the 
non-relativistic v1 = jfj^, which vanishes because the bulk modulus vanishes. In general 
we can use the fact that the energy density depends only on the density and obtain = 

dP\ — f dU\~^ dP _ K . 
dU\s ~ \dp I dp ~ pdpU- 



pdpln (dp(p^/^^{p) 



+ 1.76 



»(!■ 



(4.16) 



The density-independent result in the large-density limit is dictated by causality. It would 
be very interesting to compare these results with the ones obtained from the normal modes 
in the gravity side as in [8, 9]. 



4.3 Subleading Terms 

To study the subleading terms, we obtain the difference between the low-temperature limit 
for F from equation (4.4) and the full numerical result that we can obtain from the methods 
in [12]. 

To do so, we proceed as follows: We fix 'I'o and then vary p. Beyond the dominant 
constant term, there is some small temperature dependence in p that appears in the nu- 
merical results at OiT^). This is compensated for by using the actual numerical value 
/3(^'o, T) in (4.5) and compare this to the numeric result for the same \I'o and T. From that 
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Figure 5: The coefficients in tfie subleading term in tlie free energy parametrized as SF = —aT'^ — 
bf^ (left) and 6F = -af^ (rigfit). On tlie left, the dots indicate a regime in which the fits for the 
smaller term are not numerically reliable anymore. 



point on, we ignore the small variation in p and consider it fixed. It turns out to be most 
convenient to consider a temperature range over which the significant subleading term in 
F is suppressed by at least 10~^ and p varies at most by a factor of 10^^. Because the 
numerical accuracy in F is approximately 10"'', the systematic errors from ignoring the 
variation in p appear then roughly at the same level as the numerical noise. 

It turns out that this subleading term is always negative and interpolates from ap- 
proximately 5F — ^^j, i.e. ~ ^5-, at small density to precisely 5F = — at large 
density. The latter is consistent with the results in the massless case in [12] and [8], from 
which we expect at large densities cy = 2 ^ . The most convenient ways to parametrize the 
numeric results are 5F = —oT^ — bT^ and 6F = —aT^, which are not distinguishable at 
the level of the straightforwardly achievable numerical accuracy and are shown in figure 5. 
An exponential suppression of the subleading term, that would be suggestive of an energy 
gap, is however ruled out. 



5. General cases 



Even though it is not the focus of this paper, let us briefly comment about what happens 
in other cases than the 2+1 dimensional defect. In principle, the methods discussed in sec- 
tions 3 and 4 are straightforwardly generalized to fundamental matter supported in (d+1) 
dimensions. In the following (•)^ and (•)^ will denote dimensionless quantities analogous to 
the definitions in sections 2 and 3. The setups that come most straightforwardly to mind 
are the usual supersymmetric D3-D7 and D3-D3 setups that are T-dual to the D3-D5 case. 



12 3 

t X y z 
D3 X X X X 

D7 X X X X 



or 



background 
probe (1 + 1) 
probe (3 + 1) 

with the induced metric and the action given analogous to equation (2.2) and (2.4), re- 
spectively. For the action in the form of equation (2.8), we have to keep in mind that in 



5 6 7 8 9 



(5.1) 



X 

XXX 
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Figure 6: Left: The boundary value ^'o as a function of the density p for various setups. The large- 
density approximations are shown dashed. Right: $^ as a function of the density prf, illustrating 
also the p ^ 1 limit. 

a D3-Dp embedding, the factor (1 — 4'^) is replaced by (1 — xi/2^(p-fi-i)/2^ which is in our 
supersymmetric cases (1 — ^f^)*^/^. The solutions for the gauge field are 

duAdtiu) = = (5.2) 

/p2u2d + (l_^(-^)2)(p-d-l) 



but we will not review further details of the backgrounds as they are readily found in the 
literature [10, 19, 9] and they are not essential here. The most important property that 
we need to know for the the embeddings ^ at this point is that near it — >• 0, in the D3-D7 
system ^ ~ m-iU + qu'^ and ^ ~ fhiu + ciu In n in the D3-D3 system. 

To obtain the brane profiles, the radial coordinate has to be defined by = ^7^, and 

the large-^rf expansions of are ^0 - y^^o(1 - *o)^ ^^"^ ^0 - ^^o(l - ^o)- We show 
the corresponding ^'o as a function of the density pd = d^J^\^^^Q in figure 6. 

We find that at large densities, ^0 scales proportional to p^^^'^- As in the 2+1 case, 
there seems to be an exact relation C3 = and ci = pi^o- 

To compute the thermodynamics, we first note that the overall factor p^^"^ in the free 
energy in equation (4.4) becomes In analogy to (4.4), the factor in the integral 

coming from the AdS structure is / and in the large-^ limit, the solution (5.2) gives us 

a term in the Legendre-transformed action such that we can reproduce the ^ £]y/p,co[ 
part of the integral (4.4). Hence, the free energy becomes 

= pf^'^^'MPd) - Pdf + 0(p(i"'^)/'^) (5.3) 

and the entropy is just the baryon number as in the 2+1 dimensional case. 

5.1 3+1 Dimensions 

To discuss ^diPd) let us start with the 3+1 dimensional system. In this case, the counter 
terms from [23] are structurally the same as in the D3-D5 case, and hence the integral is 
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similar to (4.6), up to the appropriate powers in the expression. The result is shown in 
figure 6 in comparison to the 2+1 defect. 

We find that the power scaling at small densities approaches ^sips) ~ ^373- This 

means that the leading behavior in p3 is = const. — P3T3, i.e. the free energy (and also 
the internal energy) per unit volume (given by the meson mass scale) is proportional to 
the mass - independent of the baryon number density. Essentially, this implies that the 
configuration is unstable in this regime, with a negative pressure proportional to — M^g^^^. 
As the chemical potential is monotonously increasing (i.e. we have a positive "density of 
states") this instability is different from the one studied in [24]. Furthermore, there the 
"funnel" on the gravity side is becoming very narrow, whereas here it has narrow but finite 
width and is becoming very long. Still, ^' is becoming very small and hence the tension 
of the brane that falls into the horizon is close to the tension of an appropriate number of 

fundamental strings attached to a brane that does not extend down to the horizon. 

- _— 1/3 

The subleading term in F3 is just pg , giving again just the baryon number times 
the meson mass scale as we observed in the 2+1 system. 

At large densities, $3 becomes constant, close to 1, so the free energy density is ap- 
proximately F ~ p{p^^^ — T) which is consistent with the picture of an induced length scale 
discussed in section 4.2 and mimics again a Fermi surface. The sound speed in this limit 
is now Cy = \ as in [25, 26]. 

5.2 1 + 1 Dimensions 

The D3-D3 system is profoundly different. While the overall prefactor of the integral in 
(4.6) and first counter term in (4.2) just follow dimensional arguments, there are a few 
logarithmic divergences. Hence, the appropriate counter term from [23] for the scalar is 
now ^^/7^(e)^(l + l/ln(e)). In contrast to the 3+1 and 2+1 dimensional systems, where 
the boundary term gave us the appropriate variation of the action w.r.t. the scalar, c5rh, 
the variation is now —rh6c. Hence, we have to do a "Legendre transformation" and add 
the extra boundary term cm to the action in order to consider a fixed mass. This counter 
term cures the divergence coming from the scalar but there remains one divergence coming 
from the logarithmically divergent At ~ /i + lne. While the variation of the action is finite, 
and also the Legendre transformation can be trivially made finite by using p{At — pine), 
the action itself is divergent. 

In principle, one can add extra boundary terms to cure this divergence. For example in 
[9], the authors use a counter term ^^^^2X0 e ^ which has, however, a non-trivial variation 
^ ^Ay. In the limit e — )• and for a finite variation 5^^, this counter term is 

identical to the Hawking- Ross boundary term [27], and cancels precisely the boundary term 
in the variation of the gauge field. This is because, just as with the scalar, implicitly this 
term also includes a Legendre transformation, and in terms of the electric flux (or density) 
the asymptotic variation of the gauge field is bAt ~ p In e, so it diverges at the boundary and 
results in a finite variation pbp. Coincidentally, this counter term is identical to adding 

~2 

^ In e and then doing the finite Legendre transformation - or first doing a (diverging) 

-2 

Legendre transformation pAt and then adding — ^ Ine. 
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For practical aspects, the integral for <I> does converge only logarithmically in the UV 
limit, so we have to take care of this using approximate analytical solutions to extend the 
numerical result to even smaller x. We show the result for $1 in figure 6. There, we see that 
$ converges to approximately 0.596 at large densities pi ^ 1, giving us physics with an 
induced length scale similar to the other cases, and now a speed of sound of = 1, which 
is again the "conformal" result. At small densities or large masses /Oi <C 1, approaches 
$1 ~ which implies again that the energy of the system is described by a pressureless 
gas of mesons. Similarly, the subleading term giving rise to the pressure implies that that 
the pressure is at least of order p^, but numerical accuracy implies that we cannot give a 
more detailed result and also cannot exclude a very small but finite (positive or negative) 
pressure at vanishing density. 

One might worry in how far the redefinition u pi = ^1 in the logarithmic terms causes 
some mixing between the density and the temperature. As the final action however is free 
of logarithmic divergences, there should be no such remaining terms. Furthermore, we see 
that our results for the free energy rely certainly heavily on using the correct boundary 
terms. Further discussions of peculiarities of this 1+1 dimensional system can be found in 
[9], which focuses on this case, albeit not in this limit. 

6. Conclusions 

In this paper, we used holography to explore the low-temperature regime of fundamental 
matter with finite mass and baryon number density coupled to the usual Nc ^ 1, A/" = 4 
SYM theory above the deconfinement phase transition. The fact that we used a top- 
down approach guarantees us that the field theory is consistent and well defined, giving 
more relevance to our findings. Using Nf Dp probe-branes in a background of Nc ^ Nf 
D3 branes, we found a new temperature-independent scaling solution of the probe brane 
embeddings. On the field theory side, this manifests in the fact that the low-temperature 
physics is only governed by the density to mass ratio, p/M^ and interpolates between a 
mass-dominated low-density regime and a density-dominated large-density regime. On the 
gravity side the radius-dependent size of the compact sphere of the probe brane geometry 
forms a long "funnel" of finite size, with the scaling solution in the asymptotic "opening". 

First, we concentrate on a defect supported on 2+1 dimensions that is dual to a D3-D5 
setup. From the usual thermodynamic relations, we obtained the properties of this matter, 
which suggest that there is a new kind of quantum liquid. The entropy is just given 
by the baryon number times Nc\^/2 or the quark number times which implies a 

large ground state degeneracy. Consistent with this, the chemical potential has a negative 
term —TNcy/\/2, that is consistent with a Bose-Einstein distribution with a degeneracy 
proportional to the baryon number. 

The heat capacity however vanishes completely in the scaling solution, which means 
that it is of higher order in temperature than a classical Fermi or Bose gas. To explore 
the possibility of a gap, we explored the subleading terms in the free energy numerically. 
It turns out that the heat capacity is not exponentially suppressed ~ so there is no energy 
gap - but rather interpolates between cy oc at small densities and cy oc at large 
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densities; in the latter case the coefficient matches with the analytic results of the massless 
case in [8] and [5]. 

In the small-density limit, the internal energy is just given by the baryon number times 
the baryon mass scale, or precisely the overall quark number times the quark mass. At the 
same time, the pressure, bulk modulus and sound speed vanish; suggesting that the system 
becomes a pressureless gas of quarks. This is, however, still not a classical non-interacting 
gas, as the leading terms in the pressure are of higher order in temperature and pressure 
than for a classical gas, i.e. they vanish more quickly as T — )• and p — )• 0. 

At large densities, the internal energy is given by the quark number times an energy 
scale that is proportional to yj~p with a corresponding leading term in the chemical potential 
proportional to y^. This is, however, not to be interpreted as a Fermi surface, but rather as 

a length scale of the order of \j ^ that is induced by the finite density at strong coupling 
and determines the ground state energy. Note that this is much shorter than just the 
geometrical separation of quarks. Consistent with this, the pressure and bulk modulus are 
proportional to N^p^j^-^ and the sound speed becomes times the speed of light. The 
cross-over between the high and low density regimes occurrs when the induced length scale 
is of the order of the quark mass. 

We also looked at the cases of fundamental matter supported in 3+1 (D3-D7) and 1+1 
(D3-D3) dimensions. The embeddings and the general structure of the thermodynamics 
are very similar, but there are also some important differences. In the 3+1 case, the 
system exhibits a new thermodynamic instability characterized by a negative pressure at 
small density, and in the 1+1 case it turns out that there are some non-trivial logarithmic 
divergences. Applying our analysis to this case and adding additional counter terms, such 
as the one discussed in [9], yields qualitatively similar results to the 2+1 dimensional case. 

To study the properties of this quantum liquid a bit further ~ in particular in the 
2+1 defect case - it would be interesting to obtain properties that are obtained from 
linear response theory and two point functions, such as the viscosity, thermal and electric 
conductivity, quasiparticle spectrum and also the speed of sound to confirm consistency of 
the results. On the gravity side, these are straightforwardly obtained from the equations 
of motion for perturbations of the world-volume gravitational and gauge fields. Certainly, 
there may be some mixing between the different modes, but dealing with coupled systems 
in probe-brane backgrounds both analytically and numerically is straightforward [5, 9, 28]. 
Furthermore, the equations should simplify by following the limiting procedure considered 
in section 3. 

In the 3+1 dimensional case, it would be interesting to further explore this new insta- 
bility, and it would be also interesting in both the 3+1 and 1+1 cases to study the physics 
and subleading terms more thoroughly and to study dimensional peculiarities in the 1+1 
case. 

Certainly, it would also be interesting to see in how far some of the properties that we 
find are realized in experimentally studied systems. 
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